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Abstract 

A closed system of the equations for the local Bloch vectors and spin correlation functions is 
obtained by decomplexification of the Liouville-von Neumann equation for three magnetic qubits 
with the exchange interaction, that takes place in an arbitrary time-dependent external magnetic 
field. The numerical comparative analysis of entanglement is carried out depending on initial 
conditions and the magnetic field modulation. The present study may be useful for analysis of 
interference experiments and in the field of quantum computing. 
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INTRODUCTION 



The phenomenon of entanglement (inseparability) is one of the properties of quantum 
systems. A quantum entanglement is the key source of such applications as superdense 
coding, teleportation, quantum calculations researched by the quantum information theory 
0, 0,y]. The most known and natural way creating the entanglement is the global unitary 
evolution caused by the interaction between subsystems, which is controllable by external 
fields. The aim of this work has been to make a comparative analysis of the measures of 
entanglement in a three-qubit system located in a magnetic field, because up to now contrary 
to two-qubit system, there has been no unified measure entanglement for three qubits. 



THE MODEL HAMILTONIAN 



The Hamiltonian of three coupled by exchange interaction magnetic qubits e,p,n (par- 
ticles with spin 1/2) placed in an external variable magnetic field H = (Hi, H 2 , H 3 ) looks 
like 

H = m + W + W + 2 J ep s^ + 2J e X*r + 2J p "«, (1) 

where h\, h p , h™ are the Cartesian components of the external magnetic field multiplied by 
gyromagnetic ratio of the corresponding qubits; = \oi ® cxo ® <7o, s p = ~<7o <8> <Ji ® oq, 
s" = |o"o <8> cr ® &i is the matrix representation of spin operators of magnetic qubits; the 

Pauli matrices are equal ^ ocr o = ^g^^) cr i = ^^Q^) cr 2=^^ (/) ,CT3= (o ^1 

® is the symbol of direct product j3|; J ep , J en , J pn are the constants of isotropic exchange 
interaction between qubits; the summation over e,p,n is absent. 



DECOMPLEXIFICATION OF THE LIOUVILLE-VON NEUMANN EQUATION 

The Liouville-von Neumann equation for the density matrix p, describing the dynamics 
of a three-qubit system, looks like 

id t p=[H,p], p(t = 0)=p . (2) 
Let us present the solution of the equation (2) as 

P = Q-fiVyfa ®crp® cr 7 , p + = p, R 000 = 1, Spp = 1. (3) 
o 

Hereinafter summation is taken over the repeating Greek indices from zero up to four, and 
over Latin indices from one up to three. The three coherence vectors (the Bloch vectors) 
widely used in the magnetic resonance theory, are written as 

Rm = Sppai ®a ®a , (4a) 



Rom = Sppa ®ai®a , 



Rooi = Sppao ®a ®ai. 



(4b) 
(4c) 
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They characterize the local properties of individual qubits, whereas the following tensors 

Rkqo = Sppa k <g) a q ® o-Q, (5a) 



Rkoq = Sppa k ® cro ® cr, 



Rokq = Sppa ® cr fc <g> cr, 



(5b) 
(5c) 



Rkqi = Sppa k <g> a 9 <g> CT; (5d) 

describe the spin correlations. 

The length of the generalized Bloch vector b is conserved under unitary evolution. 



b = ( 6 ) 

The Liouville-von Neumann equation accepts the real form in terms of the functions R a p-y 
as closed system of 63 differential equations for the local Bloch vectors and spin correlation 
functions 

dtRqOO = ZilqhiRlOO + £mlq(J eP RlmO + J^RlOm), (7a) 

£-mlq 

{J ep R mm + J pn R olm ), (7b) 

dtRoOq = Silqh^Rooi + £lmq{J en Rlom + ^ Rolm) , (7c) 

dtRqkO — ^Uqh^RlkO + £ imkh P RqmO + J £P ^kmq{RmOO ~ Romo) 

+ J en £lmqRmkl. + J W> '^ImkRqmh (7d) 

dfRqOk — ^ilqh^RlOk + Zimkh^RqOm + J en £qmk(RoOm ~ RmOo) 

~^~J P ^lmqRmlk H~ ^ImkRqlmi (7g) 

dtRoqk = £ilqh P Roik + £ imkhf Roqrn + J Pn £qmk(RoOm ~ Romo) 

~^~J P ^ImqRlmk J ^ImkRlqmi (7f) 



dtRqkl ^imqh^Rmkl ^imkh^Rqml ^irnl^i Rqkm J P £-kmq(RmOl Roml) 

+ J en 6 qm i(Rokm — Rmko) + J Pn £kml(RqOm — Rqmo) (7g) 

for the set of initial conditions. 

In system (7), assuming, for example, J en = 0, J pn = 0, we get the closed system of 
equations for the description of two-qubit dynamics 

dtRqo = ^Uqh^Rio + e m i q J ep Ri m , (8a) 



dtRoq — Sii q h p Roi + e m i q J ep R, 
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(8b) 



dtRqk — ^ilqh^Rlk + Eimkh^Rqm. + J eP ^kmq(RmO ~ Rom) 5 (8 C ) 

where R q0 = Sppa q (g) cr , -Ro? = SjWo <8> cr g , -Rfc<? = SjWfc <8> cr 9 . 

The concrete calculations will be carried out for the following initial conditions. 

The fully separable state (S) is 



the biseparable state (BS) is (if 



S>=\111>, (9) 



BS >= 001 > + | 010 >), (10) 



V2 

the Greenberger-Horne-Zeilinger maximally entangled state (GHZ) is 



V2 

the Werner entangled state (W) is 



GHZ >= —=(\ 000 > + | 111 >), (11) 



| W >= -={\ 001 > + | 010 > + | 100 >), (12) 
v3 

the mixed state (Mix) is 

-i 

Po = x\ GHZ >< GHZ | +— ^(1 W >< W \ + \ V >< V \), (13) 

where 

| V >= -^(| 110 > + | 101 > + | 011 >), 1/3 < x < 1. (14) 
v3 

The length of the generalized Bloch vector (6) for pure states (9 -12) is equal to y/7, and for 
the mixed state (13) at x — | it is equal to v3- 



MEASURES OF GLOBAL ENTANGLEMENT IN THE THREE-QUBIT SYSTEM 

We shall define the measures of entanglement of three qubits 0] on solutions of the system 
(7), by introducing the tensors of two-partcle entanglement: 

m ijO — RijO ~ RmRojOi (15a) 
m>iOj = Rioj — RmRooj, (15b) 

The tensors (15) are equal to zero, if the two-particle correlation functions (5a, 5b, 5c) are 
factorized in terms of the Bloch vectors (4). With the help of these tensors we shall define 
the measure of two-particle entanglement in two-qubit system as 

m = rriij, (16) 
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where mij = Rij — R^R^j. The three-particle entanglement tensor is obtained by 
subtraction from three-particle spin correlation function R^f. (5d)of all products of the Bloch 
vectors and the tensors of lower order 

m ijk — Rijk — Ri00 m 0jk — RojO m iOk — RoOk m ijO ~ RmRojoRoOk- (17) 

This decomposition can be considered as a kind of cluster expansion in statistical mechanics 
In terms of these tensors the measure of three-particle entanglement becomes Q 

m SM = m% k . (18) 

This measure is equal to zero when the tensor R^k can be expressed through two-particle 
entanglements (15) and the local Bloch vectors. It is also applicable for pure as well as for 
mixed states. 

The concurrence |§| depends on all six reduced density matrices, and for pure states it can 
be represented as 



1/9 1 

C*3 = — ^= y 6 - [- + R m + Rl iQ + i?00i + + RmOn + ^Omn)]' (19) 

The measures tusm and C3 are not normalized to 1. 

The measure of global entanglement 0| i n the form ^lj is expressed through the reduced 
matrices of individual qubits p e = Sp pn p, p p = Sp en p, p n = Sp ep p according to the formula 

m B = l- g|o+^o+jooi (20) 

For the initial GHZ state of three qubits the three-tangle measure, expressed through the 
Cayley hyperdeterminant is equal to 

m K = 4pnp 88 , (21) 
where the diagonal density matrix elements determine population probabilities 

Pn = o (-^soo + -R030 + -R003 + -R330 + -R303 + -R033 + -R333 + 1) ; (22) 

o 

P88 = o (--^300 — -^030 — -^003 + -^330 + -^303 + -^033 — -^333 + 1)- (23) 
o 

The measure of global entanglement is not equal to zero for fully inseparable pure states 
[l3| and is determined by the expression 



m L = y(l- i4o)(l " 0(1 " Rm)- (24) 

This measure is equal to one for initial GHZ state and to 8/9 for the W state, and it takes 
the same values for the measure. 
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NUMERICAL RESULTS 



Let us consider two cases, when the external magnetic fields dependent on the dimen- 
sionless time r = ut: 

the resonant ( relative to e qubit) circularly polarized field (R) is 

H = — ( COST, SITIT, — ), (25) 

UJ UJ OJ 

and the nonresonant circularly polarized field (NR) is 

H = — ( COST, SlUT, — ), (2d) 

UJ UJ UJ 

where uj is the frequency of the external field , — and ^ are the dimensionless amplitudes of 
transverse and longitudinal field correspondingly. Let us assume, that the fields, operating 
on the qubits e, p and n, are equal to h e = H, h p = 2H and h n = 4H, respectively. 
To perform the numerical simulation, we have chosen the following parameters in units of 
2tt x 100 MHz, uj — ujq — 1. It corresponds to the longitudinal field H 3 = 2.3487 Ts for the 
proton resonance of e qubit. The exchange constants are equal to J ep = —0.2; J en = —0.1; 
jpn _ _q 3 anc [ (J j 1 = q_3 j n same units. In figures 1 and 2 the solid line corresponds to 
the NR field, and dashed line - to the R field. 

The S state. The measure tusm is not sensitive to the NR field (amplitude of fluctuations 
is less than 0.001), but it shows a resonant behavior in the R field (Fig.l). The measures 
C3,rriB,mL display irregular fluctuations with a low amplitude in the NR field and with a 
high amplitude in the R field. Fig. 2 shows the dynamics of measure m^. It can be seen 
that the behavior of measures msM and is similarly. 

The BS state. The measure tusm has a periodic character in the NR field and shows 
irregular fluctuations with a high amplitude in the R field. The amplitude of fluctuations in 
the R field is 3 times larger than in the NR field. The entanglement measures C^^mB^vn^ 
display irregular fluctuations with a high amplitude, and they are also little different for the 
R and NR fields. The measure C3 is shown in Fig. 2. 

The GHZ state. The measure msu is well distinct for R and NR fields. The measures 
C3, vtlb, irijCi TTT'L exhibit low-amplitude fluctuations in the NR field. In the R field the en- 
tanglement is less than in the NR field at all times, but the occurring fluctuations have a 
high amplitude. Fig. 2 shows the dynamics of measure tuk- 

The W state. The measure tusm has an oscillatory character at times of < r < 15 and 
also it is also practically indiscernible for R and NR fields. For r > 20 the amplitude of 
fluctuations in the R field is 3 times higher than in the NR field (it is not shown in Fig. 1). 
The entanglement measures Cs,mB,mL display a qualitatively similar behavior in both the 
R and NR field. Fig. 2 shows the dynamics of measure mg. 

The Mix state. It can be seen from Fig. 1, that the measure msM insignificantly changes 
in comparison with the initial value in the NR field, while the R field appreciably reduces 
the entanglement. 

As a matter of fact, the solution of equation (2) for a constant external field looks like 

p = exp(—iHt)p exp(iHt). (27) 

Evidently, if [po, H] = 0, then p = po for all t. Thus the value of an arbitrary function 
dependent on solutions of equation (2) will be equal to the initial value. In other words, the 
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Figure 1: Entanglement dynamics msM versus initial conditions 



density matrix will be fixed for this Hamiltonian. The influence of time-dependent magnetic 
field is insignificant for uj » uji, therefore all entanglement measures are constant and 
equal to the initial entanglement. 
In the model considered, a decrease in the transverse field amplitude leads to the decrease 
in the amplitudes of fluctuations for all entanglement measures, regardless of the initial 
conditions. 

We shall consider the influence of qubits e,p as two fluctuators 3| o n the n qubit 



dynamics. From system (7c) it follows for J en = 0, J pn = 0, that the spin flip probability of 
the free qubit n from an initial state -Rooi(O) = 0, -Roo2(0) = 0, i? 03(0) = 1 is equal to [131 

P = 1^91. (28) 

Fig. 3 shows the spin flip probability of n qubit in the resonant field and the dependence of 
this probability on the field deformation caused by the presence of qubits e,p in the initial 
S state. It can be seen that the Rabi oscillations transform to beats. This is in qualitative 
accordance with the results of work 3 • 
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Figure 2: Entanglement dynamics of the measures mi, C3, mj{, mg versus initial conditions 



CONCLUSION 

The closed system of equations for the local Bloch vectors and spin correlation func- 
tions for three magnetic qubits with the exchange interaction, that takes place in any time- 
dependent external magnetic field has been derived. The numerical comparative analysis of 
entanglement measures has been made, depending on initial conditions and the magnetic 
field modulation. 

From the experimental point of view, the measure vhk is more acceptable as it is ex- 
pressed through the populations pn, p 88 . The measures C 3 and m SM for GHZ conditions 
which also show sharp difference for R and NR fields but they depend in a complicated 
manner on all elements of the density matrix p . 

The present study will permit one to bring closer the theoretical results 0,0, El 11, HI, HI 
to the possibility of experimental corroboration. 

The proposed approach can be realized without any specific difficulties for the four-qubit 
system the detailed dynamics of which is described by 255 equations with the length of the 
generalized Bloch vector for pure states b = ylE. 
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Figure 3: The influence of e, p qubits as fluctuators (a thick line) on spin flip probability of the 
free n qubit n (a thin line) in the R field 

The author is grateful to Zippa Anna Anikeyevna for fruitful discussion and constant 
invaluable support. 
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